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Abstract 

The Stieltjes constants 7fc(a) appear in the coefficients in the regular part 
of the Laurent expansion of the Hurwitz zeta function C{s,o,) about its only 
pole at s = 1. We generalize a technique of Addison for the Euler constant 
7 = 7o(l) to show its application to finding series representations for these 
constants. Other generalizations of representations of 7 are given. 
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Introduction and statement of results 

The Riemann zeta function has but one simple pole, at s = 1 in the complex plane 
P,[l9l[22]. In the Laurent series about that point, 

c(^) = ^+f:^P^(^-l)^ (1-1) 

S — 1 ,n K\ 

7fc are called the Stieltjes constants [ll[5l[IIl[T3l[l6l[I71[T8l[2Tl[23], where 70 = 7, 
the Euler constant. These constants have many uses in analytic number theory, 
asymptotic analysis, and elsewhere. Among other applications, estimates for 7„ may 
be used to determine a zero-free region of the zeta function near the real axis in the 
critical strip < Re s < 1. The Stieltjes constants can be related to sums over the 
complex zeta function zeros, as well as to the Li/Keiper constants (e.g., [7]). 

The Hurwitz zeta function ({s, a) may also be analytically continued to the whole 
complex plane C/{ 1 } . It has the Laurent expansion 

r^f \ ^ , ^ (-l)^7fc(a) ^ ^ 

C(s, a) = — y + }^ [s - 1) , s^l. (1.2) 

fc=0 

Again C{s,a) has a simple pole at s = 1 with residue 1. By convention, one takes 
7fc(l) = 7fc [U]. In this paper, we are concerned with certain series representations of 
the constants 7fc(a). 

We let Pi{x) = Bi{x — [x]) = x — [x] — 1/2 be the first periodized Bernoulli 
polynomial, and {x} = x — [x] be the fractional part of x. Addison [2j gave an 
interesting series representation for 7 

1 T 00 2"-l T 1 00 2" -I 

7 = -+-yn y — - — =i-Lyn y -—^ — -. (1.3) 

2 2„t'i ^i^..2m(m + l)(2m + l) 2^, ^^^_,^^m{2m - 1) ^ ^ 



His approach uses an integral representation for the zeta function in terms of Pi. We 
note that various transformations of Addison's result are possible [31 [9l [20]. In this 
paper, we describe several generalizations of Addison's result. We indicate how a 
family of summation representations for 7^ may be obtained. Representative in this 
regard are the first two Propositions and their proofs. 

We let i/j = r'/r be the digamma function [Ij, and = J2p=i^/p the nth 
harmonic number. We have Hn = ip{n + 1) + with 7 = 
Proposition 1. We have (a) 
2^-1 r(^-_i) 



i=l m=2i- 



2 m(m + l)(2m+l) 



Inm ^ ln(m + 1/2) ln(m + 1) 



m 



2m +1 



m + 1 



[lA) 



and (b) 



3^-1 



-7i 



j=l m=33- 



(j-l)ln3 



m(3m + 1) (m + l)(3m + 2) 



Inm ^ 31n(m + 1/3) ^ 31n(m + 2/3) ln(m + 1) 



m 3m + 1 3m + 2 m + 1 

Proposition 2. Let A; > 2 be an integer. Then we have 



:i.5) 



-7i 



00 fe^ — 1 

Ei E 

j=l m=ki- 



2km + A; + 1 
2km{m + 1) 



+ il){km) — il){km + k) 



1 



(j 



In k 



k 

-^E 



=1 L 



k 



1 + In (m + ^) 1 + In (m + 



km + 



1 



km + 



:i.6) 



Alternatively, on the right side we may write the difference 



fc-i 



ip{km) — ip{km + fc) = — 



1 



mk + p 
3 



km—l 



km+k—l- 



[1.7) 



Series representations that generalize Addison's approach may also be obtained 
for 7fc(a). As an illustration, we have 
Proposition 3. Let Re a > 0. We have (a) 



OO 1 oo 



1 1 °° 1 
7n(a) = — In aH 1 — — 



2a 4 ^0 4- ^t^o {a + hj){a + h + bj){2a + b + 2bj) 



and (b) 



oo oo 



In a In^ a I ^-^ „ ^-^ 

7iw = — -^+-E2"^E 



2a 2 4 



u=0 j=0 



ln(a + bj) ^ ln(a + b + bj) ^ln[a + b{j + 1/2)] 



a + bj 



a + b + bj 



2a + b + 2bj 



b=2- 



(1.9) 



Part (a) is a known result [23]. However, we obtain it in a different way. 
From Proposition 3(b) we obtain 
Corollary 1. We have 



7i 



=27ln2 + ln2 2-7i. 



:i.io) 



From the proof of Proposition 3(a) we obtain 
Corollary 2. An integral representation corresponding to the Addison series (1.3) is 
given by 



1 , 1 /-i/l-xX^ a'^.i. 



2 2Jo \l + x 



:i.ii) 



n=l 



Our next Proposition generalizes Lemma 2.5 of [3]. We have 
Proposition 4. For any integer n > 1, we have (a) 



7 +7i 



n 



VI — x" 1 — X 



E 

fe=i 



A; In n + In In 



X 



x"" -'dx, (1.12) 



(b) and for Re a > 0, we have 



In a — ipi^a) = In a + 70(a) 



7lna+-ln^a-7V'(a)+7i(a) = ~ [j 



n 



1 — X 



1 \ 00 



k=l 



In n + In ( In ( — 

X 



(1.14) 



TT 



TT 



7^ + — In a + 7 In^ a + - In^ a - 7^ + — V^(a) + 2771(a) + 72(a) 



6 



6 



n 



1 



VI — x" 1 — X 



E 



k=l 



k"^ In^ n + 2/c In n In ( In 



In' In 



3^ diOc • 
(1.15) 

Our last Proposition generalizes relation (2.87) of [6] for 7. We introduce the 
constants ([6j, Proposition 11) 



1 /•! 

Pn+l = r / {-X)ndx 

nl Jo 



n 



(1.16) 



.=1 ^+1 

where (z)„ = r(2; + n)/r{z) is the Pochhammer symbol and s{k,£) is the Stirling 
number of the first kind. Then we have 
Proposition 5. Let Re a > 0. Then we have (a) 



In a — ip{a) = In a + 70(a) = (n — 1)! 



Pn+l 



n=l 



and (b) 



n~l 



-ln2a + 7i(a) = EP"+iE(-1)' 

^ n=l fc=0 



- l\ ln(fc + a) 
y J {k + a) 



(1.17) 



(1.18) 



From part (a) of this Proposition we obtain the following, where 2^1 is the Gauss 
hypergeometric function. We have 



5 



Corollary 3. For Re a > we have the representation 

lna + 7o(a) = -/ — r- — r . (1-19) 

a Jo [In^ (1^) +n^] V 

Recent analytic results for the Stieltjes constants may be found in [5] and [6], and 
a convenient and computationally effective asymptotic form in jl4j. The work [6] 
includes an addition formula for the constants together with series representations. 
Although the Stieltjes constants have been investigated for approximately a century, 
many open questions concerning them remain, especially including characterizing 
their arithmetic properties. 

We recall some relations useful in the following. 

From the representation (cf. e.g., [22], p. 14), 

1 1 r°° 

C(s) = + --S x-'-'+^^Pi{x)dx, Re s > 1, (1.20) 

s — 1 2 Ji 

we obtain 

= + (-1)"" r ^ " - '-^'"^ r ^ " 

(1.21) 



From (1.1) we have 



lim 



-l)"7n. (1.22) 



(s - l)«+i 

Equation (1.12) extends to 

a^^ a^"'* /■°° Pi(x) 

((s,a) = — + --S ' M+i ^^. Res>-1. 1.23 

2 s — 1 Jo [x + ay+^ 



We have the connection to differences of logarithmic sums 



7i(a) -7i(&) = 

n=0 L 

where a,b e C/{-l, -2, . . .}. 



In-* (n + a) In-' {n + b) 
n + a n + b 



J>1, 



(1.24) 



Proof of Propositions 

Proposition 1. (a) We have from (1.12) and (1.13) the integral representation (cf. 
12], p. 5) forn > 1 



In 



[ Pi(x) ^ — (n — logx)dx. 
Ji x^ 



(2.1) 



We let f{x) = -Pi{x) and put [2] 



92{x) = f{x) - -f{2x), 



(2.2) 



so that g2 is the rectangular function 



g^ix) = 1/4, < {x} < 1/2 



-1/4, 1/2 < {x} < 1. 



(2.3) 



We have from (2.2) 



E 

n=0 



g2{2^x) 



E 

n=0 



/(2"x) f{2''+^x) 



2' 



n+1 



Therefore, from (2.1) we obtain the representation 



In 



°o 1 ^OO 

E- 



(?2(2'^X) 



2"^ Ji x^ 



^ x{n — In x)dx 



(2.4) 



7 



^ln"-i(2--|/)[n-ln(2-'^l/)]ci|/, 



(2.5) 



with the change of variable y — T'x. We write the integral over y as dy — 



Y^JLv S2j" then interchange the resulting double sum in (2.5). Therefore, we 



obtain 



oo J „2^' 



7n 



EE 

j=0 i/=0 



92{y) 



2 J y 



^ ln"-i(2-^|/)[n-ln(2-^|/)]dy. 



We now speciahze to n = 1: 

oo J r2^^'g,{y) 



-EE/ 

oo j „ 

= -ee/ 



[1 - ln(2-'^y)]dy 



[1 + i/ln2 — lny](ij/ 



i=0 "^2^ 



'4>[l + ^ln2-l.,]d,. 



2J 2/ 



We then use the integral 



Jm y^ 4 \Jm Jv 



.+1/2 J 



dy 



We find 



oo 



1 + lnm 4(1 + ln(m + 1/2)) 1 + ln(m + 1) 



m 



2m + 1 



+ 



m + 1 



(2.6) 



(2.7) 



(2.8) 



-7i = TEi E 



j=l m=2i- 



l + i2^1n2 



1 + lnm 4(1 + ln(m + 1/2)) l + ln(m + l)' 



m(m + l)(2m + l) m 



+ 



2m + 1 



m + 1 
(2.9) 



With the use of partial fractions on the first term on the right side of this equation, 
we obtain the form (1.4). 
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For part (b), we proceed similarly, employing the function 



so that E^=o 



g3(3"a:) 



gsix) = fix) - -fi3x) = 1/3, 0<{x}< 1/3, 



= 0, 1/3 < {x} < 2/3, 
= -1/3, 2/3 < {x} < 1, 
= f{x). We use the integrals 

Jm y'^ 3 \Jm im+2/3y 



(2.10) 



and 



/ 

Jrr 



3[m(3m + l) (m + l)(3m + 2) J ' 



(2.11a) 



1+2/3 / y^ 



1 + lnm 3(1 + ln(m+ 1/3) 3(1 + ln(m + 2/3) ^ 1 + ln(m + 1) 



m 



3m + 1 



3m + 2 



m + 1 



(2.116) 

Again, an application of partial fractions gives the final form (1.5). 

Proposition 2. Proposition 1 is further extended with the use of functions for 
integers k >2 

9k{x)^f{x)-^f{kx), (2.12) 

with Z^^o ^''^k^^^ — /(^)- Since / is periodic, Qk is also periodic of period 1. Prom 
the expression 

9k{x) = ^ (l - ^) + li^^} - {^}' (2-13) 



9 



we obtain the values of gk(x) on subintervals ior j — 1,2, . . . ,k. We have 

\j - 1 i 



A; 'A; 



(2.14) 



This enables the development of summation expressions for any 7„ as a result of 
integrations over a function gk{x). We have from (2.1) 



1 r gkik'x) 



u=0 



1 



In" x{n — lnx)dx 



-E r^ln"-^(fc-^y)[n-ln(fc-'^y)]d|/ 



^ln--\k-'y)[n-ln{k-'y)]dy. 



j=0 ;/=0 



So at n = 1 we obtain 



y 



-71 = EE 

j=0 ;/=0 



^fc(y) 



[l-Hk-'y)]dy 



''"9k{yl^l + ll^k-\ny]dy. 



ly y^ 



i=o -"^ 

Now using the piecewise values in (2.14), we have 



/ 

Jrt 



m y 

k 



k 



k 



m+e/k dy 

m+{e-i)/k 



E 

£=1 



k 



k 



{km + £){km + i- 1) 



2km + k + l 
2km{m + 1) 



+ ip{k'm) — ip{k'm + k). 



We also have 



/„ 



1 - 



1\ (^-1) 
k 



L 



m+l/k dy 

ln|/— 
m+{e-i)/k y^ 



(2.15) 



(2.16) 



(2.17) 
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k 

E 

e=i 



1 - 



k 



-1) 
k 



l + ln(m + ^) l + ln(m + f) 



km + £-l 



km + i 



combining (2.16), (2.17), and (2.18), we arrive at the Proposition. 
Proposition 3. (a) From (1.13) and (1.2) we have 



hm 



C{s,a) - 



S-1 



— Ina + 



2a 



-/ 

^0 



°° Pi (x) 



{x + a} 



;dx, 



so that 



fix) 
{x + a)' 



1 f 
7o(a) = -lna + — + 

Using the functions gk of (2.12) we have 

Jo {x + a)2 ^ ~ ^0 k" Jo 

~ E l.2v E /. 
u=o'^ 3=0-''' 

We now speciahze to k — 2 and obtain 
/ 



- 1 ^ fffc(y) 
.=0- j=oJi ik-'^y + a)^ 



X + a)2 
dy. 



dx 



a; + a)2 .-^ 4^^ ^ 



^■+^ 92iy) 



-1 OO -| oo 



•j+1/2 «+l 



, {2'-y + a)2 



'j+i/2y (2-y + a)2 



oo T oo 



1 



4 S S (« + ^i) (« + ^' + ^'J) (2a + 6 + 26j) 



6=2- 



1 oo 

-Y 



2^ 



2a + 6 
26 , 



a + 6 



6=2- 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



(2.22) 



We now apply the functional equation and duplication formula for the digamma 
function: 



JO 



! A ^ 



ix + a)2 



v=0 



2^ 



'2a + 6' 
2h . 



-2^ T -- 



6=2- 



11 



-I oo 

-y 



y=0 



6=2-" 



6=2-" 



1 

2a 

1 

2a' 



(2.23) 



We employ a standard integral representation for ip ([lOj, p. 943) to write 

,1 j^^2a2"-l _ ^a2"-l^l 



E 

!/=0 



6=2- 



T 

lim 

i/=0 



t- 1 



- In 2 



lim 



1 ^^a2'^+i-l _ 



■rft- (T + l)ln2 



lim 



t - 1 

7/^(a2^+^) - V^(a) - (T + 1) In 2 
= — '?/'(a) -|- In a, 



(2.24) 



where the above interchange of summation and integration is justified by absolute 
convergence of the integral, and where in the last step we applied the large argument 
asymptotic form for ip (|1], p. 259). Therefore, from (2.20), (2.23), and (2.24), we 
find (1.6). 

(b) Using (1.13) and (1.2) we have 



lim 



C'(s, a) — a^' 



In a 



givmg 



7i(a) 



(s-l)2 



In a In a 



In^a In a Pi{x 



2 2a Jo (x + a 



— \l — \n{x + a)]dx, 
I 

(2.25) 



+ ^ + 



Piix) 



2 ' 2a ' Jo (x + a)^ 
With k = 2m (2.12) we have 

'■^ Plix) r. , . , ^^ , ^ 1 >^ 92{y) 



[1 — ln(x -|- a)]dx. 



(2.26) 



r -™-[l - ln(x + a)]dx = V - y / 



v=0 ^ j=o • 
12 



J (2--y + a)- 



■[l-\n{2--y + a)]dy 



-1 OO 1 oo 

-T-T 

^ i/=0 ^ j=0 



[1 - ln(2-"|/ + a)\—-^^-—.. 



(2.27) 



We then carry out the integrations and simphfy. Using (2.25) completes part (b). 
Corollary 1. We put a = 1/2 in (1.9) and manipulate the resulting series to write 

11 2 



7i(l/2)+ln2+iln^ = i^f:2-+^f: 



i/=0 



+ 



l + 2-''+ij l + 2-^+i(j + l) l + 2-^+i(j + l/2) 



-1 oo oo 



ln(l + bj) ^ ln[l + b{j + 1)] _ ^n[l + 6(j + 1/2)] 



1 + bj 



1 + 6(j + 1) 



l + 6(j + l/2) 



6=2- 



-1 oo 



ln(l + 2j) ^ ln[l + 2(j + l)] _ ^Infl + 2(j + 1/2)] 



(2.28) 



l + 2j ■ l + 2(j + l) l + 2(j + l/2) 

By the method of (2.23)-(2.24), the first line on the right side of this equation has 

the value — In 2(7 + In 2 — 1). Prom Proposition 3(b) at a = 1, the second line on the 

right side of (2.28) is 71. The following evaluates the third line of the right side of 

(2.28). 

Lemma 1. We have 



-1 00 



ln(l + 2j) ^ ln[l + 2(j + l)] _ ^Inil + 2(j + 1/2)] 



l + 2j 
Proof. We write 



l + 2(j + l) 



l + 2(j + l/2) 



1 






~ 2 


71 





+ W2. 
(2.29) 



-1 00 

4 ^ 



ln[2(j + 1/2)] ^ ln[2(j + 3/2)] _ 2ln[2(j + 1)] 



j + 1/2 



J + 3/2 



i + 1 



ln2^ 



j=0 



+ 



1 

+ 4 



71 I 2 



7i + 7i ( 2 



7i 



j + 1/2 i + 3/2 j + 1 
(2.30) 



where we used (1-13) at j — 1. Then the sum of (2.30) is given by 



!^(21„2-l) + i 



7i o + 7i 



271 



13 



In2, , , 1 

2 n2 - 1 + - 

2 ^ ' 2 



71 ( 2 ) ^ 71 + In 2 



71 



71 



+ In' 2. 



(2.31) 



In obtaining this result we have used 7^(0 + 1) = 7fc(a) — (In^ a) /a, that follows easily 
from the series form oi C,{s,a)^ k = 1. 

Then combining the terms on the left and right sides of (2.28) yields (1.10). 

Corollary 2. We have, by applying for instance Theorem 2.2 of [3], 



1 /I 



Vl + a; 



x^dx 



p + 



J + Hp/2 - iJ(p_i)/2, Re p > -1. 



(2.32) 



Then, with the interchange of the summation and integration justified by the absolute 
convergence of the integral, we have 



1 /I -x 



V 1 + X 



J2x^''-'dx = -l + J^[H, 



p/2 



(p-l)/2jp=2"-l5 



n=l 



n=l 



(2.33) 



giving 



1 1 /I - X 

2^2/0 \TTx 



E 

n=l 



X 



2"-l 



00 

dx = -y^ 
2 ^ 



n=l 



P+l 



p=2"-l 



(2.34) 



By the method of (2.22)-(2.24), the Corollary follows. 
Proposition 4- Part (a) is based upon the integral 

This equation follows easily by term-by-term integration, as we have 



' (-^ + In' (- ] dx 

\inx 1 — xj \x, 



mx j^Q-'^ 



14 



OO -| 

= -r(.) + r(. + i)X:^-^^, (2.36) 

where we used the change of variable u = — Inx. We apply Lemma 2.4 of [3] that we 
write in the following form. For any integer n > 1 and Re a; > 0, we have 



lux 1 — X ^ ra^ ^j=o 
We then proceed as in the proof of Lemma 2.5 of [3J to find 



Vlnw \ — u) \u) Jo Vl — x" 1 — X 



X 



OO 



fc=l 

(2.38) 



We perform ^ on (2.35) and (2.38), and (1.12) follows. 
Part (b) is based upon the integral 

(-^ + x^-Hn' (-] dx = T(s)[sC(s+l,a)-a-'], Res>-1, Re a > 0. 

Jo Vina; 1 — xJ \xj 

(2.39) 

Using Lemma 2.4 of ^ we find 

/ V— + u--' In^ (-]du= n —] In^ ( -) y n^^x-'^"-^dx. 

Jo \mu l — uj \uj Jo \1 — x" 1 — X/ Vx/ ^ 

(2.40) 

Taking s = in (2.38) and (2.39) gives (1.13). Performing (^)"' , j = 1,2 on 
(2.39) and (2.40), respectively gives (1.14) and (1.15). 

Proposition 5. We employ the generating function ([6j, (2.78)) 
OO II 

Y: Pn+iz"-' = - + — r, \z\ < 1, (2.41) 

^1 z ln(l - z) 



15 



so that by (2.39) we have 

r(s)[sC(s+ l,a) - a""] = Y.Pn+1 C x^'-^l - xf-Hn' (-\ dx. (2.42) 

n=l •'^ ^-^^ 

For part (a), we take s = in this equation and evaluate the integral using the Beta 
function. For part (b), we differentiate both sides of (2.42) with respect to s and put 
s = 0. We find 

-7ln a- ^In^ + 77/^(0) -7i (a) = f^Pn+i x^-^l-xY-^ln {\n dx. (2.43) 

By logarithmic differentiation we obtain the integral for Re s > — 1 

r >n' (i) 1. (in (i)) = + 1) - + 1)1. (2.44) 



riving at s = 



'x'l„flnfi)V.= -Il±M^. (2.45) 



/o \ \xj J + 1 

Then by binomial expansion we have 

"v-', , ./n-l\[7 + ln(fc + a)] 



^.-(l-.)-ln(l„(-))d. = -i;(-l) 



/c / k + a 



fc^ \ k J k + a 
By using (2.43) and (1.17) we find (1.18). 

For Corollary 3 we apply the integral representation (2.85) of [6] for Pn+i, so that 
by (1.17) 

(n - 1)! /■°° du 



lna + 7o(a) = E^^VT^/ 
~, {a)„ Jo 



^1 {a)n Jo (l + M)"(ln'M + 7r2) 
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1 Z"*^ 1 1 / 1 \ 

= - / ^ ^7 T 2F1 1, 1; 1 + a; du. (2.47) 

a Jo (In^M + vr^) (1 + m) ' H l + u) ^ ^ 

Remarks. From (2.6) we have for n > 1 

-7n = f:E r"^E7''7^V-^1^2)^l^"~'~'^[^ + '^1^2-lny]rfy, (2.48) 

and more generally from (2.15) for A; > 2 we have 
^ f'''^' 9khj) (n - 1' 

With (2.20) and (1.8) we have recovered a known integral representation for InF 
or the polygamma functions ip^^'^ in terms of Pi (e.g., [8J, Section 6.3). However, our 
route has been very different from say the use of Euler-Maclaurin summation. 

The result (1.10) also follows from the relation (^(s, 1/2) = (2* — l)(^(s). The 
method of Corollary 1 shows an instance of the functional relation contained in (1.9). 

In Proposition 5(a), when a = 1, we have (1)„ = ra!, and we recover (2.87) of [B]. 
Other forms of part (b) may be written by using the integral representation 

du 



/n-l\ ln(fc + a) --"^ 
to ^ \ k ) k + a 



Tl — 1 ) ! 

—— — '- 2Pi(l, 1 - n; 1 + a; u) 

a L, a 



Inw 
(2.50) 



This expression has been obtained by inserting an integral representation for In on the 
left side, interchanging summation and integration, and making a change of variable. 

We note that representations of other constants may be obtained from Lemma 2.4 
of [3]. We let li(z) = Jq dt/ Int be the logarithmic integral. Then we have 
Corollary 4. For any integer > 1 we have 

n 1 X 00 



x" ^ 

k=l 



-Hl-y) + Hy) = l [^^--^jE^^' -'dx. (2.51) 
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Our results for 7^ also imply representations for the coefficients rij that appear 
in the Laurent expansion of ('/( about s = 1 (e.g., ^). As an example, we have 
771 = 7^ + 271. 

Summary 

We have shown that the method of Addison may be generalized in several direc- 
tions. It applies not only to the Euler constant 7, but to the Stieltjes constants 7jt for 
the Riemann zeta function. It further extends to 7^(0) for the Hurwitz zeta function. 
Moreover, one may develop a set of parameterized series representations by using 
the step-wise functions gk{x). On yet a broader scale, the series developments are 
applicable any time the first periodic Bernoulli polynomial Pi is integrated, including 
in numerous occurrences of Euler-Maclaurin summation. 

With Proposition 4, we have generalized representations of 7 given by Ramanujan 
and Berndt and Bowman [3]. 
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